On collapse in the nonlinear Schrodinger equation with time dependent nonlinearity. 

Application to Bose-Einstein condensates. 
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It is proven that periodically varying and sign definite nonlinearity in a general case does not 
prevent collapse in two- and three-dimensional nonlinear Schrodinger equations: at any oscillation 
frequency of the nonlinearity blowing up solutions exist. Contrary to the results known for a sign 
alternating nonlinearity, increase of the frequency of oscillations accelerates collapse. The effect 
is discussed from the viewpoint of scaling arguments. For the three-dimensional case a sufficient 
condition for existence of collapse is rigorously established. The results are discussed in the context 
of the meanfield theory of Bose-Einstein condensates with time dependent scattering length. 

PACS numbers: 03.75.Lm, 03.75.Kk, 03.75.-b 
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Stimulated by achievements in experimental real- 
ization and management of Bose-Einstein condensates 
(BECs) , considerable attention has recently been paid to 
possibilities of stabilization of condensates by using time 
variation of the scattering length |2jj3j, Li • By analogy 
with the famous Kapitza pendulum [5(, which acquires 
an additional dynamically stable point when the pivot 
is rapidly oscillating, it has been suggested 0, 0] that 
nonlinearity varying rapidly enough and changing sign 
during oscillations can stabilize a quasi-two-dimensional 
(2D) condensate. According to 1] stabilization can be 
achieved even when inter-atomic interactions are charac- 
terized by a negative definite scattering length. 

Previous studies were based either on qualitative ar- 
guments, like variational approach 0, El using the 
Gaussian ansatz and more accurate moment method , 
or on direct numerical simulations of a multidimensional 
nonlinear Schrodinger (NLS) equation. While most of 
the obtained results report similar conclusions about the 
critical collapse, approximate character of the exploited 
approaches results in some discrepancies in conclusions 
about 3D collapse. In particular, for the negative mean 
(i.e. averaged over one period) scattering length the au- 
thors of Ref. were not able to arrest the collapse, while 
stable solutions were reported in Ref. 0]. 

Thus, the present situation of the theory clearly 
demonstrates lack of exact results. It turns out that rig- 
orous statements, which constitute the main goal of the 
present Letter, are available in a case of a sign definite 
scattering length 01 . Being subject to a number of con- 
straints our results do not solve the problem completely, 
but allow one to understand the effect of the scattering 
length modulation on the solutions of the 2D and 3D NLS 
equations. In particular, we prove that variation of neg- 
ative definite scattering length with any frequency gen- 
erally speaking does not arrest collapse, i.e. one can al- 
ways choose an initial condition blowing up at finite time. 
Qualitatively, the above statement can be conjectured on 
the basis of earlier results. Indeed, as it has been proven 
in 0, dissipation cannot arrest the overcritical collapse, 



but only changes the sufficient condition of its existence. 
Meanwhile it is known that the varying nonlinearity in 
the NLS equation can be transformed into a time depen- 
dent dissipative term (see e.g. |7j). Thus one can expect 
that, in a general situation, time dependent nonlinearity, 
even rapidly varying, will not arrest collapse, but only 
change conditions for this phenomenon. In this Letter, 
we formulate sufficient conditions for the collapse in the 
case of the negative mean scattering length and find that 
in the 3D case oscillations of the nonlinearity are favor- 
able for the collapse, in the sense that increase of the 
frequency of oscillations leads to decrease of the upper 
bound for time of collapse. 

In the case of a BEC in the mean-field approxima- 
tion, the problem is described by the Gross-Pitaevskii 
(GP) equation 0], which in the absence of the external 
trap potential is also known as the NLS equation, for 
which the problem of collapse was intensively studied for 
a long time Jsj. Taking into account that the parabolic 
trap potential (the typical one for the most experimen- 
tal settings) does not affect the fact of the existence of 
collapse phenomenon in the case of a constant scatter- 
ing length JlfJl , we restrict our considerations to the NLS 
equation. Besides reproducing the main qualitative re- 
sults for the BEC, such a statement provides generality of 
the results, as the NLS equation is the well known model 
equation for numerous physical phenomena. In particu- 
lar, the results described below for the critical case are 
directly applicable to the problem of beam focusing in a 
stratified Kerr medium. 

Statement of the problem and scaling arguments - Let 
us consider the dimensionless NLS equation 



dt 



= -VV«-fl(w*)|^«|Va 



(1) 



where tf) u = ip u (x,t) and x € M. D with D = 2,3 be- 
ing the spatial dimension. The nonlinearity coefficient is 
considered to be varying with a period T = 27r/w, i.e. 
g(t) = g(t + 2w) and to be bounded and positive definite: 
.92 > g{t) > gi > for all t. It will be convenient to 
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introduce the notation ip(x,t) = tpi(x,t) for the solution 
of Eq. with u> — 1 which in this way does not contain 
any free parameters. Thus, ip u (x, t) — ^/uj(p (uit, \/ujx) is 
a solution of (JIJ for a given u, whenever <p(t, x) solves Q 
with u> = 1. 
The energy 



- ff (wi)|^ 



If, (2) 



and the number of particles N u — J [ip^dx play a spe- 
cial role in the analysis of the blow up phenomenon (if not 
specified, hereafter the integrals are taken over ~R D ) @. 
One easily verifies the following relations 



The last equation, as well as the link between the so- 
lutions ip and ip u , mean that existence of a blowing up 
solutions of Eq. with u> = 1 and negative energy 
implies existence of a blowing up solution of Eq. at 
any oscillation frequency of the nonlinear term. If that 
happens in the critical case (D = 2) the collapse occurs 
with the same number of atoms, while in the 3D case the 
number of particles required for the collapse decays as 
1/a/o7 as the frequency goes to infinity. 

When the nonlinear term is a positive constant (i.e. 
when g(t) =const> 0), the solution of Q blows up at a 
finite time, provided the energy is negative Below we 
will show that this also happens in the case of varying 
nonlinearity, where energy will be required to be nega- 
tive at the initial moment of time. As in Eq. Q the 
link between initial energies of the solutions ip u and (p is 
given by E u (0) = cj( 4 ~ d >' 2 Eo (hereafter we simplify the 
notation introducing Eq — E u (0)) which means that by 
increasing the frequency one increases the modulus of the 
energy of the blowing up solution ^^(x, t) proportionally 
to ui and y/uj in the 2D and 3D cases respectively. 

In the case at hand, however, the energy J5J) is not a 
constant any more but is governed by the equation 



dK, 
dt 



(4) 



The energy grows during half periods with dg/dt < 0, 
and thus in principle may acquire positive values, even 
being initially negative. Thus the rigorous results of the 
NTS collapse cannot be applied straightforwardly. They 
can however be modified to provide sufficient condition 
for the collapse of solutions of Eq. 0), which we will dis- 
cuss in the next two paragraphs (here we concentrate on 
physical applications of the theory; more mathematical 
details and generalizations will be published elsewhere). 

"Early-time" collapse - Let us start with the most sim- 
ple, but allowing rather general considerations, situation 
where g{u>f) is growing during the first half-period, and 
thus the energy is decaying. Then, E u (t) < for the 
interval t € [0, T/2] (provided the solution exists) and to 



get a sufficient condition for the collapse it is enough to 
require that it happens during the first half-period (that 
is why we call it "early-time" collapse). 

This can be done by a slight modification of the stan- 
dard arguments [|J . To this end we introduce the quanti- 
ties Y(t) = J \x\ 2 \i/}w\ 2 dx, and Z(t) = Im J x-Vip^ ^dx, 
which solve the equations: 



dY(t) 

dt 
dZ{t) 
dt 



= -4Z(t), 

= -DE u (t) + (D - 2) J |VV>a 



dx 



(5) 
(6) 



From © and © it follows that, if dg/dt > 0, Eq < 0, and 
Z > (hereafter Y(0) = Y , Z(0) = Z Q ) one can obtain 
the estimate Y(t) < 2DE t 2 - 4Z t + Y , from which it 
follows that the blow up occurs at a finite time < Tq < 

Z 



oo, where To 
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2DE Q 



Imposing now 



the condition To < T/2 we obtain a requirement for Yq: 

Y a <Y* = D\E \T 2 /2 + 2Z T. (7) 

This condition and the requirements Eq < and Zq > 
constitute the sufficient conditions for the collapse to 
happen during the first half-period. 

The obtained result has transparent physical mean- 
ing. Indeed, compared to the standard, time indepen- 
dent problem, a new condition J7J appeared. Since Y(t) 
is a mean squared width of the wave packet, the new 
condition requires the initial wave-packet to be localized 
sufficiently well to decrease the blowing up time, making 
it less than the first half-period. 

Combining the above result with the scaling arguments 
of the preceding paragraph, one concludes that for any 
oscillation frequency of the nonlinearity with initially 
positive derivative, one can find an initial condition for 
collapse at finite time, in both 2D and 3D cases. 

Sufficient conditions for the collapse in the 3D case 
- Condition J7J looses its practical sense in the case of 
rapidly varying nonlinearity, i.e. when T — > 0. Then, in 
physically relevant situations, collapse cannot occur dur- 
ing the first half-period and one has to consider a more 
general situation which will be restricted to the 3D case. 
Since the sign of the energy is of primary importance and 
assuming that initially the energy is negative, Eq < 0, in 
order to establish a sufficient condition for the collapse 
we have to control the change of the energy E UJ (f) in time. 
We will do that using the ideas due to Tsutsumi 0. 

Taking into account that g{ojf) is a periodic function 
with a period T, we consider an interval t € [T„_i,T„], 
where T n — nT with n being an integer, and assume that 
the solution exists in this interval (more precisely in the 
interval t 6 [0, T n ]). As we have shown in the preceding 
paragraph, the way how the nonlinearity is changing dur- 
ing the first half-period is relevant for the early-time col- 
lapse. Now we relax this constrain, and choose the most 
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"unfavorable" for collapse (because of initial grows of the 
energy) situation, where dg/dt is time definite on each of 
the half-periods, with dg/dt < for t € {T n _ x ,T n - T/2) 
and dg/dt > for t e (T n - T/2, T n ). 
Next we define two functionals 



£{t) 



lwg 2 -7sWkU 4 )^ 



(8) 



Integrating by parts e~ at E u (t) with respect to time and 
using Q we obtain for t > t\ and for some positive con- 
stant a > the following relation 

e- at E u (t) = e^E^ih) -a f e- as £(s)ds . (10) 

Let now £ > <i and t, h € [T„_i,T„ - T/2). Then one 
has £ (ti) < E(ti) and Eq. IjlOl) allows us to obtain 

d 

/ e Q(t - s) £(s)ds < e^-' 1 '^!) . 
The last formula implies 



e- as £(s)ds<0 ifT(ii)<0. 



(11) 



ti 



We have assumed that initially the energy is negative, 
i.e. i?o < 0. Then, using the continuity arguments, which 
take into account that ljll|l is valid for all t and t\ from 
the interval [T n _i, T„ — T/2), we obtain that in the first 
half period £{t) < 0. Next we observe that E u (T/2) = 
£(T/2) < and that E^{t) is a decreasing function in 
the second half period. Hence E^iT) < 0. Noting that 
£{T n ) = E u (T n ) and £(T n - T/2) = T W (T„ - T/2) for 
all n for which the solution exists and applying the pre- 
vious arguments for the first n periods, we deduce that 
the initial condition To < guarantees that T W (T„) < 0. 
In other words periodically varying nonlinearity with def- 
inite sign cannot result in a change of the sign of an 
initially negative energy. 

For the next consideration we recall JSJ and JSJ, rewrit- 
ing the last expression for the 3T> case as follows : 
dZ/dt — —2£(t) > — 2£ (t) where £ (t) is a continuous 
function defined by: £{t) — £{t) when t € [T„_i,T n — 
T/2] and £{t) = E when t € [T n - T/2,T„]. Then the 
following estimate for Y(t) holds 



Y{t) < Y + 4 



ds 



Z + 2 £(a)da 



(12) 



Let us define r = r(t) through the relation t = nT + r, 
where n is chosen to be the largest integer assuring that 
nT < t and thus < r < 1. Then, the first integral in 



(|1H is trivially computed, while for the second one we 
obtain 



t r s rnT 

£(a)dads< (nT — s)£(s)ds + \Eq\ T < 



^0 







E T<[?--^)+\E \T 



The last formula and Eq. I|12() allow us to obtain the 
estimate as follows 

Y(t) < 2E T 2 n 2 - (5E Q T 2 + 4Z Q T)n + Y Q + \E \T(13) 

From this inequality we can find the number n» deter- 
mining the latest period during which blow up occurs 
(at that number the right hand side of the inequality be- 
comes negative). In this way we obtain that the blow up 
occurs at t < T„ r , where 
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Thus we have outlined the proof of the following 

THEOREM: Let ip u be a sufficiently smooth solution 
of 0) in 3D case, the initial condition for which is char- 
acterized by Eq < and Zq > 0; then blow up occurs at 
a finite time t < T», where T* is given by 

It is worth to emphasize that although we considered a 
situation where change of the scattering length is initial- 
ized with the "negative" half-period of dg/dt, the above 
estimates obviously applies for any initial value of dg/dt. 

Estimates for real condensates - Let us now discuss 
the qualitative picture emerging from the obtained re- 
sults, restricting the analysis to the 3D case. We notice 
that the temporal characteristics of the collapse are rel- 
evant to the theory of a BEC with periodically varying 
scattering length due to experimental constraints on the 
frequency of the oscillation of the nonlinearity, emerging 
from the fact that change of inter-atomic interactions in 
practice is achieved by means of the Feshbach resonance, 
controlled by varying external magnetic field. The same 
physical phenomenon can result in creation of molecules 
from pairing atoms, in originating excited atomic states, 
etc. The respective processes are not described by the 
meanfield GP equation (the NLS equation), what re- 
stricts the range of meaningful frequencies. On the other 
hand, relevant frequencies are bounded from below by 
characteristic times of the condensate's life. 

Although the sufficient condition gives only an upper 
bound for the time of the collapse, we will treat the quan- 
tities T) and T» as the collapse times (conjecturing that 
in a general situation decrease/increase of each of these 
quantities results in decrease/increase of the time of the 
collapse). Then the first observation is that Tq < T*. 
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Second, the upper bound for collapse T* decreases as the 
frequency of oscillation grows. This is in sharp contrast 
to what is predicted in the 2D case with the sign alter- 
nating scattering length 0. 0- E| The third feature to 
be mentioned is the dependence of the collapse time not 
only on the number of atoms involved but on the aspect 
of the initial distribution. 

To connect our results with realistic experiments we 
provide the estimates using the data from Ref. [i"l"|| . 
where observation of the collapse of a BEC controlled 
by Fesbach resonance (with monotonically changed mag- 
netic field) was reported. We consider a cloud of con- 
densed 85 Rb atoms initially having a Gaussian distribu- 
tion normalized to the number of particles N and char- 
acterized by the radius r (in dimensionless variables): 

H x ,°) = rvl'^n ex P (~f^)> which S ives z o = 0, 

Y = Wr 2 /2, and E = - 2 , /2 % r3 . For the en- 
ergy to be negative in the described situation, one must 
have N > N cr — 2 1 / 2 6ir 3 / 2 r. 

Change of the scattering length is modeled by the for- 
mula a s (t) = a°(l — B uy_ Bo ), where a° = — 20.1 nm, 
the position of the resonance peak is Bq = 154.9 G and 
the width of the resonance is A = 11 G. We consider the 
initial magnetic field B(Q) — 166 G, which corresponds 
to the initial scattering length a s (0) = — 0.18nm, and 
the amplitude of the field oscillations 10 G (which for the 
frequency 1000 Hz corresponds to the speed 6.37 G/ms of 
the change of the magnetic field). Considering the ini- 
tial radius of spherically symmetric cloud to be 16.5 /jm 
(which corresponds to r = 1 in the dimensionless units) 
one obtains that the link between iV and the real num- 
ber of particles AT is given by M (N/7) ■ 10 4 (the 
unit of the dimensionless time corresponds to 0.116 s), 
and thus Af should exceed Af cr = 67498. In the case 
at hand the "early collapse" happens at (physical) times 
bounded by to ~ 13.77/vW — JV cr s for the frequencies 
lu 2 < to 2 w 0.756 • 10- 2 (W - N cr ) Hz. If frequency in- 
creases, or the scattering length initially decreases, col- 
lapse occurs at later times, bounded by T*. Although the 
respective analytical expression for T* is readily obtained, 
it appears to be more informative to present dependence 
of the upper bound of the collapse time vs the frequency 
of the scattering length graphically (see Fig. ^| . 

Conclusion - It has been established that nonlinearity 
periodic in time but sign definite does not prevent col- 
lapse in two- and three dimensional condensates with a 
negative mean scattering length. A sufficient condition 
for the collapse has been formulated which implies the 
possibility to create initial configurations of a condensate 
which will blow up in a finite time. 

The sufficient condition of the 3D collapse is obviously 
not the optimal estimate for the time of the collapse. 
This is not only due to the fact that in the course of the 
proof some estimates were shortened (and lower precision 
was the price to pay for that), but mainly because the 




FIG. 1: Upper bound for the collapse time, vs frequency for 
different numbers of atoms. The frequency, loo, defined in the 
text for each of the curves is 4.3492 Hz (solid line), 13.0428 Hz 
(dashed line), and 17.9253 Hz (dotted line). 



proof does not involve a specific law of the variation of 
the nonlinearity (but only the fact that it is periodic). 
The respective improvement of the estimate, as well as 
its generalization to the critical case (considered in the 
present Letter only for the case of early collapse), are left 
as open questions. 

To conclude, the above results can be directly gener- 
alized to the nonlinear Shrodinger equations with higher 
nonlinearity and to the NLS equation with a periodically 
varying dissipative term. 
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